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We study the flavor dependence of the medium modification factor Rpp for pp collisions for scenario 
with formation of a small-size quark-gluon plasma (QGP) for RHIC {^/s = 0.2 TeV) and LHC 
{^/s = 2.76 TeV) energies. We And that at pr ~ 10 GeV the pion spectrum is suppressed by 
~ 20 — 30(25 — 35)% for RHIC (LHC), for D (B) mesons the suppression effect is smaller by a 
factor of ~ 0.7 — 0.8 (0.5). The flavor hierarchy Rpp < Rpp < Rpp is held at pr < 20 GeV for RHIC 
and at pr < 70 GeV for LHC. This gives a sizeable reduction of the heavy-to-light ratios of the 
nuclear modification factors Raa as compared to that in the standard scenario without the QGP 
production in pp collisions. 

PACS numbers: 


I. INTRODUCTION 

There is general consensus that the observed in ex¬ 
periments on AA collisions at RHIC and LHC suppres¬ 
sion of the high-pT particle spectra (jet quenching) is 
due to radiative and collisional parton energy loss in 
the QCP formed in the initial stage of AA collisions. 
In the scenario with formation of a mini fireball of the 
QCP in pp collisions a similar mechanism should sup¬ 
press particle spectra as compared to the predictions of 
the standard pQCD that neglect the final state inter¬ 
action (FSI) effects for fast partons produced in hard 
reactions. In this case the measured pp inclusive cross 
section will differ from that predicted by the pQCD, 
dcfpert{pp —t hX)/dpTdy, by the medium modification 
factor Rpp 

_ dajpp -» hX)/dpTdy 
dapertipp -t hX)/dpTdy ’ 

where the numerator is the real pp inclusive cross sec¬ 
tion accounting for the FSI effects in the mini-QCP. Of 
course the Rpp is not directly measurable quantity, be¬ 
cause the spectrum with the FSI effects switched off is 
unknown. Nevertheless, it is clearly of great interest for 
understanding the uncertainties of the standard pQCD 
predictions due to the FSI higher twist effects. Also, 
the medium suppression in pp collisions is important for 
theoretical predictions for the nuclear modification fac¬ 
tor Raa for A A collisions, because the pp baseline spec¬ 
trum, that is necessary for evaluating the Raa, should 
account for the medium effects in the mini-QCP. In prin¬ 
ciple, the medium modification of jets in pp collisions can 
be studied experimentally by measuring the multiplicity 
dependence of the direct photon-triggered fragmentation 
functions (FFs) [I|. 

The idea that the QCP may be produced in pp colli¬ 
sions has attracted much attention in recent years (see, 
for instance, [l-Q)- It is mostly due to observation of the 
ridge effect in high multiplicity pp events at y/s = 7 TeV 
by the CMS collaboration Q, which may be caused by 


the transverse flow of the produced mini-QCP fireball. It 
is important that the conditions for the mini-QCP for¬ 
mation in pp collisions are better in jet events, because 
the multiplicity of soft off-jet particles (the so-called un¬ 
derlying events (HE)) is enhanced by a factor of ^ 2.5 Q 
as compared to minimum bias events. Therefore, even at 
RHIC energies y/s ^ 0.2 TeV the UE multiplicity may be 
high enough for the QCP formation. The scenario with 
mini-QCP in jet pp events is supported by the prelimi¬ 
nary data for pp collisions at y/s = 7 TeV from ALICE 
0 indicating that the jet FFs become softer with in¬ 
crease of the UE multiplicity. 

In we have addressed jet quenching in pp colli¬ 

sions for light hadrons within the model of jet quenching 
of 0 (see also [13-0) based on the light-cone path in¬ 
tegral (LCPI) approach 0 to induced gluon emission. 
We studied the medium modification of the y-triggered 
and inclusive FFs in ra collisions [ij and evaluated Rpp 
for charged hadrons [ll|. It was found that the medium 
effects in the mini-QCP in pp collisions may be quite 
strong, say, at pr ^ 10 CeV the spectra are suppressed 
by ^ 20 — 30% at RHIC energy -/s = 0.2 TeV and by 
- 25 - 35% at LHC energy = 2.76 TeV. 

In the present letter we study flavor dependence of the 
Rpp. It is of great interest for understanding the accu¬ 
racy of the pQCD predictions for heavy quark produc¬ 
tion in pp collisions. Also, the flavor dependence of Rpp 
is important for theoretical predictions of Raa for heavy 
flavors. The question of jet quenching for heavy flavors 
has attracted much attention in recent years (see 0 
for a short review) due to the observation of strong sup¬ 
pression of single (non-photonic) electrons from decays 
of the heavy mesons in experiments at RHIC 0 0 
(the “heavy quark puzzle”). The recent measurements 
at LHC of Raa for electrons 0 and for D mesons 0 
also show a strong suppression effect. It seems to be dif¬ 
ficult to reconcile with the expected dead cone suppres¬ 
sion of the radiative energy loss for heavy quarks pre¬ 
dicted in 01 . A subsequent reanalysis [0 of the quark 
mass dependence of induced gluon radiation within the 
LCPI approach 0 demonstrated that due to the quan- 
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turn finite-size effects, ignored in the dead cone model 
[^ . the quark mass suppression of radiative energy loss 
at low energies (< 20 — 30 GeV) turns out to be signif¬ 
icantly weaker than predicted in [2^. And at energies 
>100 GeV the radiative energy loss even rises with the 
quark mass. Calculations of Raa for electrons and D 
mesons in the scenario without the QGP formation in pp 
collisions performed in [3, [T^ using the LCPI approach 
[l^ have shown reasonable agreement with experimen¬ 
tal data. However, the experimental error bars for Raa 
for heavy flavors are very large. Also, the experimen¬ 
tal data (especially for the non-photonic electrons from 
RHIC [la, 113) restricted to relatively low px, where 
the assumption of dominance of the radiative energy loss 
and the relativistic approximation mp /Eg <C 1 for heavy 
quarks (especially for 6-quark) used in [ij, [l^ may be in¬ 
valid. Therefore, it is probably too early to conclude that 
the “heavy quark puzzle” is solved, and the flavor depen¬ 
dence of the nuclear modification factors deserves further 
investigation. In the present work we study possible ef¬ 
fect of the mini-QGP formation in pp collisions on the 
heavy-to-light ratios of the nuclear modification factors 
in A A collisions, which is connected to reduction of the 
medium suppression for heavy quarks in pp collisions. 


II. SKETCH OF THE THEORETICAL 
FRAMEWORK 


We simulate jet quenching in pp collisions within the 
approach of m- It is qualitatively similar to the scheme 
developed for AA collisions in [^, which has been suc¬ 
cessfully used in our several previous analyses of jet 
quenching in AA collisions [l3l - [lq . So only a brief outline 
of important points of our theoretical framework will be 
given here. The interested reader is directed to 0113 
for details. 

For the perturbative inclusive pp-cross section we use 
the standard formula 


do-pert (pp hX) 

dprdy 





da{pp —>■ iX) 
dplpdy 


( 2 ) 


where pt is the particle transverse momentum, y is 
rapidity (we will consider the central region y = 
0), Dfi/i is the vacuum parton—^-particle FF, and 
da{pp ^ iX)/dplpdy is the ordinary LO pQCD hard 
cross section, = pr /z is the initial parton transverse 
momentum. For the initial virtuality we take Q = plp. 
In calculating Rpp we write the real medium modified pp 
cross section, that enters the numerator of o, in a form 
similar to (I5|) 


damjpp -» hX) 

dprdy 



.,pm p. ^.daipp^iX) 

xD,/,{z,Q) 


( 3 ) 


where is now the medium-modified FF for i ^ h 

transition in the presence of the mini-QGP fireball. It 
is implicit that is averaged over the jet production 
point, and the impact parameter of pp collision. As in 
0, we used the distribution of hard processes in the 
impact parameter plane obtained using for the parton 
distribution in the transverse plane the transverse density 
distribution for quarks in the MIT bag model. 

In the standard scenario without the QGP formation 
the nuclear modification factor for AA collisions (we de¬ 
note it Raa) ™ scheme of [13] for a given impact 
parameter 6 reads 


pst (, „ „ X ^ dN{AA hX)/dpTdy 

AA[ ,PT,y) T^^^h)dap,rt{NN ^ hX)/dpTdy ’ 

( 4 ) 

where TAA{b) = f dpTA{p)TA{p — b), Ta is the nucleus 
prohle function, and the numerator of (|3]) (we omit the 
argument 6) is given by 


dGm(NN^hX) .. 

"- 1^, -- 

Here dam{NN -A- hX)/dpxdy is the medium-modified 
NN ^ hX cross section, which is given by the formula 
© for NN collisions with the medium-modified FF for 
the QGP produced in AA collision. In the scenario with 
the mini-QGP formation in A A collisions the theoretical 
Raa (which should be compared with the experimental 
Raa) can be written as 

Raa = Raa/R pp ■ ( 6 ) 

Because experimentally Raa is defined as the ratio of the 
measured AA -A hX cross section to the binary scaled 
experimental AA —>■ hX cross section, and the latter 
includes the FSI effects in the mini-QGP produced in 
A A collisions (hereafter we ignore the difference between 
the FSI effects in pp, pn and nn collisions). 

As in [I 3 , we evaluate for each fast parton tra¬ 
jectory as 

® D"), 0 Dk/dQ ), (7) 

where 0 denotes z-convolution, D^./i is the parton 
DGLAP i ^ k FF, DJJf, is the parton j ^ k FF in 
the QGP accounting for parton energy loss, and D^/j is 
FF for the parton j ^ h transition outside the QGP. 
For the stage outside the QGP for light partons we use 
for the Dh/j{Qo) the KKP 0 FFs with Qq = 2 GeV. 









3 


And for the FFs c ^ D and b 
parametrization 


B we use the Peterson 


Dh/q{z) oc 


z[l-(l/z)-eQ/(l-z)p 


( 8 ) 


with Ec = 0.06 and et = 0.006. The DGLAP FFs have 
been obtained with the help of the PYTHIA event gener¬ 
ator [ 2 ^. The medium-modified FFs were computed us- 
ing the induced gluon spectrum in the form obtained in 
1^ . We account for the (relatively small 271) collisional 
energy loss by redefining the initial QGP temperature 
in our formulas for the medium-modified FFs related to 
induced gluon emission (see [l^ for details). 

For the quasiparticle masses of light quarks and gluon 
in the QGP we use the values niq = 300 and rrig = 400 
MeV supported by the analysis of the lattice data [^ . 
for c and b quark masses we take uic = 1.2 GeV and 
nib = 4.75 GeV. We use the Debye mass p,D in the QGP 
obtained in the lattice analysis [ 2 ^, which gives ^jld/T 
slowly decreasing with T {fiu/T Ri 3 at T ^ 1.5Tc, 
/i£)/T R! 2.4 at T ^ 4Tc). 

We use (both for radiative and collisional energy loss) 
running as frozen at low momenta at some value 
For gluon emission in vacuum for this parametrization a 
reasonable choice is R 0.7— 0.8 [^[^. Since ther¬ 
mal effects can suppress the in-medium QGD coupling, 
we treat as a free parameter (which may differ for pp 
and AA collisions). It is the only parameter that controls 
the strength of the medium effects in our calculations. 

As in [T^, we use for the QGP evolution 1-l-lD 
Bjorken’s model (both for pp and AA collisions), which 
for the ideal gas model gives TqTq = T^r, where tq is the 
thermalization time. We take tq = 0.5 fm. For r < tq 
we take medium density oc r. For simplicity we neglect 
variation of Tq with the transverse coordinates. 

We fix the initial temperature of the plasma fireball in 
pp collisions from the initial entropy density determined 
via the experimental UE multiplicity density 


So = 


C dN,h 

To7ri?j dr] 


(9) 


Here C = dSjdy jdNch!dr] r 7.67 [s^ is the en¬ 
tropy/multiplicity ratio, and Rf is the typical radius of 
the created fireball. We use for Rf the parametrization of 
obtained from results of simulations of pp collisions 
performed in I™ within the IP-Glasma model. This pro¬ 
cedure gives [all = 0.2, 2.76 TeV] r [1.3,1.44] fm 

and To[^ = 0.2,2.76 TeV] R [199,217] MeV. (we 
use the ideal gas model of the QGP with Nf = 2.5). 
Note that the initial temperatures would be higher by 
~ 10 — 15% if we used the entropy from the lattice calcu¬ 
lations [al]. We ignore this fact, because in our analysis 
the temperature is an auxiliary quantity characterizing 
the QGP entropy. The entropy is this quantity that is 
only important from the point of view of parton energy 
loss (if we assume that the number density of the color 


constituents in the QGP is approximately proportional to 
its entropy density). And we determine it directly from 
the experimental UE multiplicity density. Eor AA colli¬ 
sions To was fixed using the data on the charged hadron 
multiplicity pseudorapidity density dNct/dr] from RHIG 
[ 3 ^ and LHG [s^. It gives Tq r 320 MeV for central 
Au-I-Au collisions at ^/s = 200 GeV, and Tq r 420 MeV 
for central Pb-|-Pb collisions at y/s = 2.76 TeV. 

Note that our results for Rpp are not very sensitive to 
the size of the fireball, which is not well determined. An 
analysis of the stability of Rpp to variation oi Rf per¬ 
formed in [llj has shown that ±30% change in i?/ gives 
a small effect on Rpp. This is due to a compensation 
between the rise of the energy loss with the fireball size 
(for a fixed density) and its suppression with decrease of 
the fireball density. In [ll| we argued that the stabil¬ 
ity of Rpp to variation of Rf allows to expect that the 
transveses expansion of the medium and the variation of 
the initial medium density in the transverse coordinates 
should not have a significant impact on our results for 
Rpp. 


III. NUMERICAL RESULTS 

In Fig. 1 we present the results of our calculations 
of Rpp for pions, D and B mesons for RHIG energy 
^/s = 0.2 TeV and LHC energy ^/s = 2.76 TeV obtained 
for = 0.5, 0.6 and 0.7. For pions and D mesons 
we show the results for pT > 5 GeV, and for B mesons 
for the lower limit we take px = 7 GeV to avoid possi¬ 
ble problems with applicability of the relativstic approx¬ 
imation and with treating the collisional energy loss as 
a smallperturbation to the radiative energy loss for b 
quark [l^. One can see that the suppression effect due 
to the mini-QGP formation turns out to be quite large. 
At pt 10 GeV the pion spectrum is suppressed by 
- 20 - 30(25 - 35)% for RHIG (LHC), for D mesons 
the suppression is smaller by a factor 0.7 —0.8, and for B 
mesons the effect is smaller by a factor of ^ 0.5. At RHIG 
energy all the Rpp flatten at pr > 15 — 20 GeV. This oc- 
cures due to compensation between the effects from the 
reduction with parton energy of the relative parton en¬ 
ergy loss AE/E (which increases Rpp) and from the in¬ 
crease of the effective power Ue// in the power low depen¬ 
dence of the hard cross section da/dpx oc (which 

reduces Rpp) with increase of parton energy. Fig. 1 shows 
that the medium suppression in pp collisions is not very 
sensitive to It is a consequence of relatively large 
transverse momenta kx (relative to the initial fast par- 
ton momentum) of the emitted gluons for the small-size 
QGP, because they approximately scale with the QGP 
size, Lqgp, as 1/Lqgp [HI- In [Tlj ] in calculating the 
nuclear modification factor Raa for the scenario with 
mini-QGP in pp collisions we used the value = 0.6 in 
calculation of the factor Rpp. In the present analysis we 
also use this value. 

Even though the medium suppression of the spectra 
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s'^=Q).2 TeV s'^=2.7Q TeV 



FIG. 1: Rpp of pions, D and B mesons at ^/s = 0.2 TeV 
(left panels) and ^/s = 2.76 TeV (right panels). From top to 
bottom, the curves corresponds to = 0.5, 0.6 and 0.7. 


from the QGP formation in pp collisions may be size¬ 
able, it is hardly possible to observe it directly by com¬ 
paring the experimental spectra with the theretical pre¬ 
dictions, because the uncertainties of the standard pQCD 
predictions (see for example ) are larger than the 

expected medium effect. For this reason it is difficult 
to disentangle the medium suppression and other effects 
that can change the pQCD predictions (such as modifica¬ 
tion of the PDFs and FFs, change of the renormalization 
and factorization scales) However, as was mentioned in 
Introduction, the preliminary ALICE data [l^ on the UE 
multiplicity dependence of the jet EEs for pp collisions at 
y/s = 7 TeV support the jet quenching in pp collisions. 


From Fig. 1 one sees that the flavor hierarchy < 
Rpp < Rpp is held for px < 20 GeV for RHIC and for 
Pt < 70 GeV for LHC. One can expect that in these 
PT-regions the presence of the Rpp in the formula ([6]) for 
the nuclear modification factors in AA collisions will re¬ 
duce the ratio of the nuclear modification factors Raa 
for heavy and light hadrons. However, to reach a definite 
conclusion about the effect of Rpp on the heavy-to-light 
ratios of the nuclear modification factors Raa for AA 
collisions one should account for the fact that in the sce¬ 
nario with mini-QGP formation in pp collisions the value 
of al'' should be somewhat increased to keep stable Raa 
for pions (since it is assumed to be adjusted to the ex¬ 
perimental Raa of pions). 


Comparison with the data on Raa for light hadrons 
carried out in till has shown that for the standard sce¬ 
nario the RHIC data on Raa for pions in Au-I-Au col¬ 
lisions for 0 — 5% centrality bin at y/s — 0.2 TeV sup¬ 
port ~ 0.5, and the LHC data on Raa for charged 


Au+Au, s''^=0.2 TeV Pb+Pb, s^^^=2.76 TeV 



FIG. 2: The ratios Raa/R\a (upper panels) and Raa/R\a 
(lower panels) for 0—5% central Au-|-Au collisions at ^/s = 0.2 
TeV (left) and for 0 — 5% central Pb-|-Pb collisions at yfs = 
2.76 TeV (right). The solid curves are for Raa for = 0.5 
(left) and a/A = 0.4 (right) without Rpp factors in The 
dashed curves are for Raa for = 0.6 (left) and af'’ = 0.5 
(right) with Rpp factors in ([6|. The dotted curves are for Raa 
for of’’ = 0.5 (left) and = 0.4 (right) with Rpp factors 
in ®. For the dashed and dotted curves the factors Rpp for 
light and heavy flavors are calculated with q{’' = 0.6. 


hadrons in Pb-|-Pb collisions for 0 — 5% centrality bin at 
y/i = 2.76 TeV support ~ 0.4. And for the scenario 
with mini-QGP formation in pp collisions the experimen¬ 
tal data prefer somewhat larger values: ~ 0.6 (0.5) 

for RHIC (LHC). For these values of at px > 10 GeV 
the decrease of Raa ffi® numerator of (O reasonably 
compensates the presence of Rpp < 1 in the denomina¬ 
tor of (O. This allows one to have approximately the 
same Raa in both scenarios of pp collisions. However, 
this procedure acts in a different way for heavy quarks. 
It results in some differences in the heavy-to-light ratios 
of the nuclear modification factors Raa for the scenarios 
with and without mini-QGP production in pp collisions. 

In Fig. 2 we show the ratios Raa/^aa ^aa/^aa 
for Au-fAu collisions at y/i = 0.2 TeV and Pb-fPb colli¬ 
sions at ^/s = 2.76 TeV (for 0 — 5% centrality class) ob¬ 
tained for Raa without (solid curves) and with (dashed 
curves) inclusion of the Rpp factors for 0 — 5% central AA 
collisions at RHIC and LHC. We have used the above 
values of 0.5 (0.4) for RHIC (LHC) in the standard 
scenario, and 0.6 (0.5) for RHIC (LHC) in the scenario 
with the mini-QGP formation in pp collisions. From 
Fig. 2 one sees that the mini-QGP formation can size- 
ably reduce the heavy-to-light ratios at px < 20 GeV for 
RHIC and at px < 50 GeV for LHC. At px ^ 10 GeV 
(interseting in the context of the present situation with 
the “heavy quark puzzle”) in the mini-QGP scenario the 
ratio Raa/^AA becomes smaller by ^ 20 — 30%, and 























5 


Raa/^aa by 15 — 25%. To illustrate better the role 
of the increase of q;{'’ in the mini-QGP scenario, we show 
in Fig. 2 the ratios for Raa obtained with inclusion of 
Rpp in ([5]), but for Raa evaluated without renormaiza- 
tion of (dotted curves). By comparing the results for 
two choices of one can see that the renormalization 
of Qfl’' gives a sizeable reduction of Raa/^aa’ but prac¬ 
tically does not change Raa/Raa- Thus, from Fig. 2 we 
see that the mini-QGP formation in pp collisions weakens 
the flavor dependence of the nuclear modification factors 
Raa at RHIC and LHG energies. However, the effect is 
not very strong, and due to rather large experimental er¬ 
ror bars for the available heavy flavor Raa data (Tsl - f^ 
the situation with their description in the scenario with 
the mini-QGP production in pp collisions does not differ 
significantly from our previous analysis [alii within 
the standard scenario without the QGP formation in pp 
collisions. For this reason we do not present a compar¬ 
ison with experiment in the present letter. Anyway, for 
a conclusive comparison between theory and experiment 
we need more accurate data on the Raa for heavy flavors 
at higher px- 

IV. SUMMARY 

The medium FSI effects should modify high-pr jets 
in pp collisions if the mini-QGP production occurs. We 


have studied the flavor dependence of the medium mod¬ 
ification factor Rpp in this scenario. We evaluated Rpp 
for pions, D and B mesons for RHIG (-^3 = 0.2 TeV) 
and LHC = 2.76 TeV) energies. We have observed 
that at pt ^ 10 GeV the pion spectrum is suppressed by 
- 20 - 30 (25 - 35)% for RHIC (LHC), for D {B) mesons 
the effect is smaller by a factor of ^ 0.7—0.8 (0.5) than for 
pions. The flavor hierarchy is held at 

Pt < 20 GeV for RHIC and at px < 70 GeV for LHC. We 
demonstrated that this gives a sizeable reduction of the 
heavy-to-light ratios of the nuclear modifcation factors 
Raa- At px ^ 10 GeV the ratio Raa/^aa suppressed 
by ~ 20 — 30% and Raa/^aa by ~ 15 — 25% as com¬ 
pared to that in the standard scenario without the QGP 
production in pp collisions. 
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